Probing Quantized Einstein-Rosen Waves with Massless Scalar Matter 



o 
o 

(N 



o 

o . 

I 



X 



J. Fernando Barbero gQ and Inaki GarajQ 
Instituto de Estructura de la Materia, CSIC, Serrano 123, 28006 Madrid, Spain 

Eduardo J. S. Villascnoi0 

Grupo de Modelizacion y Simulacion Numerica, Universidad Carlos HI de Madrid, 
Avda. de la Universidad 30, 28911 Leganes, Spain and 
Instituto de Estructura de la Materia, CSIC, Serrano 123, 28006 Madrid, Spain 

(Dated: April 19, 2006) 

The purpose of this paper is to discuss in detail the use of scalar matter coupled to linearly 
polarized Einstein-Rosen waves as a probe to study quantum gravity in the restricted setting pro- 
vided by this symmetry reduction of general relativity. We will obtain the relevant Hamiltonian 
and quantize it with the techniques already used for the purely gravitational case. Finally we will 
discuss the use of particle-like modes of the quantized fields to operationally explore some of the 
features of quantum gravity within this framework. Specifically we will study two-point functions, 
the Newton- Wigner propagator, and radial wave functions for one-particle states. 
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INTRODUCTION 



^ , There are reasons to believe that general relativity is not the last word as far as gravitational physics is concerned, 
(yT) ' for one thing, we still do not know how to reconcile it with quantum mechanics, the other fundamental pillar of 
] our description of physical reality. As of today we lack a complete and consistent description of quantum gravity 

■ although some proposals have made steady progress towards our understanding of the problem if not of its solution. 
' The questions to be addressed by quantum gravity must go beyond the mere finding of a consistent quantization of 

general relativity because it is possible that many of the physical concepts that work classically are inappropriate in 
the quantum setting. Let us consider, for example, the metric. Though this is the fundamental concept in classical 
general relativity, formalisms such as loop quantum gravity suggest that it may not be the basic variable at the 
quantum level. It is perfectly conceivable that the metric is some kind of semiclassical construct that emerges at 
scales much larger than the Planck scale. We are left, then, with the problem of understanding what the relevant 
vTti' object at the shortest scales is and how a metric description appears. Another concept that may not be appropriate 
rr' is that of the test particle. Its classical meaning is clear, they are point-like objects of negligible mass and, hence, 
^ with no influence as sources for the gravitational field described by a metric whose geodesies are the trajectories 
followed by these test particles in their motion through space and time. They are, in a definite sense, the tools that 
we have available to extract the non-trivial geometric content from the metric and describe spacetime physics. In the 

■ context of general relativity it is very useful and illuminating to adopt an operational point of view to define geometric 
quantities such as time intervals, lengths, and other geometrical objects with an immediate physical interpretation. 
In this way it is possible to avoid resting too much on our classical (non-relativistic) intuition and isolate the basic 
and relevant physical concepts. If a consistent quantization of gravity was available it would seem necessary to replace 
test particles by quantum objects. One of the obvious choices -though possibly not the only one- is to introduce 
quantum particles moving in the non-trivial background defined by the metric. The fact that this will in general be 
curved would manifest itself in the way wave functions spread throughout spacetime. Another possibility would be 
to couple test fields and use their quanta as test particles, yet many others may exist as we cannot be sure that the 
correct description at quantum gravity scales will be that provided by particle-like objects. 

The purpose of this paper is to explore these issues in the restricted setting provided by the linearly polarized two- 
Killing vector reductions of general relativity, the so called Einstein-Rosen waves. The general problem of introducing 
quantum test fields consistently is a non-trivial one. We will see, however, that it is possible to use some fields for 
this purpose because we will be able to solve the theory exactly both at the classical and quantum level even in their 
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presence. Specifically we will expand here the results presented in showing that Einstein-Rosen waves coupled to 
a massless scalar field can be exactly quantized. After this we will use the particle-like quanta of this scalar field as 
a way to extract information about this quantum gravity toy model in an operational way, much in the same way as 
test particles are used to derive spacetime physics in general relativity. Our concern will be the recovery of classical 
trajectories for particles and an approximate metric description. The results obtained here complement in a sense 
those already derived by considering the microcausality of this model Hi Q| that clearly show how a long distance 
limit appears in which microcausality reduces to the familiar one defined by the Minkowski metric. 

The structure of the paper is the following. After this introduction we describe in section ^ the classical solution 
of a system of linearly polarized Einstein-Rosen waves coupled to a massless, cylindrically symmetric, scalar field. We 
will use the Geroch formalism 5] and take advantage of the two Killing vector fields available in order to describe the 
model as a 1-1-1 dimensional system of coupled fields. As we will see the structure of the solutions to the field equations 
strongly suggests that the Hamiltonian of the system has a very simple form. We discuss this in section lTTll One of the 
points that merits special attention here is the treatment of the boundary terms that give rise to the Hamiltonian. In 
fact, the most salient feature of the model is the appearance of a Hamiltonian that is a non-trivial bounded function of 
the sum of the free Hamiltonians corresponding to two free, massless, axially symmetric scalars in a 2+1 dimensional 
Minkowskian background. This fact will eventually allow us to quantize the model even though it is non-trivial and 
interacting. In section llVl we apply well-known Fock space techniques, similar to the ones already used in 6], to study 
its quantization. After this we will describe several applications. We start in section by studying the two-point 
functions. They can be interpreted as approximate probability amplitudes for a particle created at a certain radial 
distance from the cylindrical symmetry axis at a certain time to be detected somewhere else at a different instant of 
time. The comparison of these amplitudes with the ones corresponding to a cylindrically symmetric massless scalar 
in a Minkowski background shows some characteristic features due to the presence of gravity such as an enhancement 
of the probability of finding the field quanta close to the axis. Microcausality is easily reviewed in this framework 
because the field commutator can be obtained from these two-point functions. As the reduced field model has two 
different fields, the "gravitational" and the "matter" scalars, we can study their commutators. We will see that they 
are not zero as a consequence of the non-trivial interacting character of the model; the commutator does not show 
the microcausal behavior that appears when only one type of field is considered. 

The main reason why the interpretation of the above two-point functions as probability amplitudes is only approx- 
imate is the fact that the states obtained by acting with the field operators on the vacuum are not orthonormal. The 
obtention of an orthonormal basis of position states in quantum field theory is a difficult problem that was essentially 
solved by Newton and Wigner for some special types of fields . By using this basis it is possible to introduce rela- 
tivistic wave functions whose modulus squared can be interpreted as probabilities to find a particle around a certain 
spacetime point. ^ We introduce analogous objects in section IVIl Due to the absence of Lorentz invariance in the 
present model some of the unpleasant features of the ordinary Newton- Wigner states are absent here. We use these 
states in a double fashion. First we study the propagator {R2\Uit2,ti)\Ri), that can be interpreted now as a proper 
probability amplitude, and compare it with the two-point function considered above. After this we define "radial 
wave functions" in terms of these Newton- Wigner states in section lVlIl and study in detail a certain family of them for 
which the evolution can be explicitly obtained in closed form as a one-dimensional integral. This provides us with the 
main tool that we have been searching for in the paper: a way to introduce quantized particles with controlable wave 
functions that we can use as quantum test particles to explore the spacetime physics of this quantum gravitational 
toy model. We will use them to study the propagation of particles in this purely quantum spacetime. 

Throughout the last sections of the paper we will make use of asymptotic techniques developed by the authors in 
to extract relevant information from the closed integral expressions that describe two-point functions, propagators, 
and wave functions. This is useful for two different reasons. In some cases it allows us to consider the "macroscopic" 
(or "classical" ) limit in which we only consider the behavior of the relevant objects at scales much larger than the 
quantum gravitational one. In others they provide us with useful tools to obtain approximate values for these object 
that are not easily derived by numerical methods. We end the paper with our conclusions and comments in section 
IVIIII and an Appendix where we discuss a useful representation for several integrals appearing in the paper. 



^ Although these so called Newton- Wigner states provide a rather complete solution to the problem of localization in some relativistic 
quantum field theories they also display some disturbing features mainly related to their behavior under Lorentz transformations 1^. 
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II. EINSTEIN-ROSEN WAVES COUPLED TO MASSLESS SCALAR MATTER 

In this section we review in some detail the classical theory of whole cylindrically symmetric spacetimes minimally 
coupled to a cylindrically symmetric massless scalar field. Although this problem has been extensively considered in 
the literature (see for instance and references therein), it is useful to gain some insight about the Hamiltonian 
formulation of the model (that we will discuss in the next section) through the analysis of the field equations and the 
Geroch formalism. As we will see the solutions to the field equations strongly suggest that the Hamiltonian of the 
system, obtained by carefully taking into account surface terms in the action principle, is a simple function of the 
Hamiltonians for two massless cylindrically symmetric scalars. We start by considering the Einstein equations for the 
system 

□'*V = o, (1) 

=8^Gjv(d0)a(d0),,. (2) 

Here R^^i^ and are, respectively, the Ricci tensor and the d'Alembert operator associated with the Levi-Civita 
connection ZJJj*' compatible with the spacetime metric g''^^^. The exterior derivative of the matter scalar field (p is 
denoted by (d^)^ and Gjv is the Newton constant. 

Whole cylindrical symmetry is characterized by the existence of a M x U{1) group of isometrics with two 
mutually orthogonal, hypersurface-orthogonal, commuting, space-like globally defined Killing vectors and a"'. It is 
assumed also that the topology of the spacetime is M"* and that the set of fixed points of the rotations generated by 
ct" = (d/da)"" defines a 2-dimensional time-like surface: the axis of symmetry {a; S K'' | ffa^'cr^cr'' = 0}. On the other 
hand the translations along the symmetry axis generated by — {d/dz)"- act freely and satisfy the elementary flatness 
condition [11.] . This condition guarantees the 27r periodicity of the axial coordinate a. With these assumptions the 
Killing coordinates z G M and a G [0, 27r) are unique up to the trivial transformations z i—t az + zq and a i-^ ±cr + (jo- 
in terms of these Killing fields the symmetry of the system reflects on the vanishing of the Lie derivatives 

Finally, we will restrict our discussion to the class of spacetimes for which the derivative {dR)a of the scalar field 
defined by 

is everywhere space-like. R is the area density of the isometry group orbits. 

In order to solve the Einstein equations (^[21) we will make use of the Geroch reduction technique 0. First, owing 
to the fact that translations have no fixed points, it is possible to rewrite the Einstein equations as equations for 
fields defined on the quotient manifold -topologically M^- comprised by the translational orbits. To do this we need 
to introduce a scalar field A and the tree dimensional spacetime metric (^^^J given by^ 

These fields, as well as </>, are well defined in the translational orbit manifold. In terms of them the equations ^ 
are equivalent to 

□<^V = -^ff<^"^'(d0)a(dlogA)6, 

□<^'A = ^5'^'«''(dA),(dA),, (3) 

Rab = ^D^a 'mb - ^(dA),(dA)fc + SttGat id^)aid^)b , 

where i?^^', Da \ and D'^' refer to 5^^^'. It is easy to prove that the rotations generated by cr" are still a symmetry of 
the reduced theory. In particular cr" is well defined in the quotient manifold and satisfies 

L'r9ab = ' L.X^O, and = . 



Here and in the following i;a = d'^t^^ a^nd ""a = S^^''' 
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At this point, although the rotations do not act freely, we can still apply the Geroch reduction to the manifold obtained 
by deleting the symmetry axis p = g^i^lcr°'<J^ — from R^, and incorporating these points at the end of the process by 
imposing regularity at the axis. In particular, outside of the axis, we can define 

P = 9ty<^'^9':y^'>Q and = ffS + V„a, . 
The equations © can then be written as 

□<=V--^.9'^""'(d0)a(dlog(Ap)), 

□ <^)A = -^g<^>'^^(dlogA),(dlog(A-V))6 
□<'V=f5'^'"'(dlogp).(dlog(A-V))6 

Ral - ^D^amb + Yp^aiAp), - ^ (d log A), (d log A)b - ^ (d log p), (d log p), + SttG^ (d(/.), (d(/.)b , 

where the notation in the previous expressions is the natural one. These 2-dimensional field equations can be solved 
in two steps. First, it is convenient to replace the field p in terms of i? = ^/\p and A to get 

□ (^)^ + g(^)afc(d^)^(dlogi?), = 

□<''logA + g<'"'^(dlogA)(dlogi?) = (4) 
□<"'i? = 

Ral - ^D^a' (di?)fc - ^ [d log A] (a [d log(A-ii?2)],) + SttGat {d4>)a{d4>), . 

Second, owing to the fact that we are dealing now with 2-dimensional field equations, D'^'i? = allows us to introduce 
a new scalar field T, the harmonic function conjugate to R, whose derivative (dr)a is everywhere time-like by means 
of 

(dT), = 6<^,'g'^'^=(di?),, 

where e'-^l is the volume element associated to g'^l- Notice that this definition is conformally invariant. It is now 
possible to use R and T as coordinates in the l-f-l reduced spacetime and introduce the fiat metric 

C = -(dr)a(dT)b + {dR)a{AR)i . 

The degrees of freedom of g^^^ are encoded in the conformal factor defined by 

(2) (2) 

9ab — ^ '^ab ■ 

In terms of these new fields, the equations I0J become 

□ + 77<^"'''(d(/.)a(dlogi?)fc = O 

□ log A + 7?<^"'^(d log \)a{A log R)b = 

r;<^>^'*[i(dlogA)c(dlogA)rf-f SvrG^ {A<P),[d^)d - (d7),(d log i?)^] r^^ 

= |a,(di?)b - 2(d7)(,(dlogi?)b) + (dlogA)a(dlogA)6 + %tiGn (d0)a(d0)fc, 

where □ is the d'Alembert operator defined by the fiat metric 77^^^. The general solution to these equations can be 
written in a very convenient form in terms of the fields 



(t>g := log A , 0s := Vl67rGAf(/> . 

By using R and T as coordinates, the Einstein equations are equivalent to two uncoupled cylindrically symmetric 
Klein-Gordon equations 



1 

'r^ 

drl - R[{dR(Pg){dT(Pg) + {dB,(ps){dT(Ps)] , 

dRl = ^[{dT(Pgf + {dR^g? + {dT<Psf + [ORCj^sf] 
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The first two are equations for a massless, axially symmetric, scalar fields in a Minkowskian background and the 
equations for 7 satisfy an integrability condition that allows us to immediately write their solution as 

7 = ^ y [(5t0<,)' + {dR<l>gf + {dTcj^sf + {dR<Psf] RdR. 
Finally the four-dimensional spacetime metric satisfying the Einstein field equations can be written as 
9ab = e^"*' [ - (dT)a(dT)b + {dRUdR)t] + ii^e"^-' {da)a{dak + e^' (dz)a(dz), . 

The form for the 'C-energy' 7 strongly suggests that the Hamiltonian of the system can be obtained from the one 
corresponding to Einstein-Rosen waves by adding the contribution of the extra scalar field. We show this in the next 
section. 



III. HAMILTONIAN FORMALISM 

In order to develop the Hamiltonian formalism, we start from the Einstein-Hilbert action in four-dimensions for 
gravity coupled to a massless cylindrically symmetric scalar (jjs- 

We have included the boundary terms needed to have a well defined variational principle. The fields are taken to 
be regular in the symmetry axis and the boundary conditions at infinity correspond to the definition of asymptotic 
fiatness introduced by Ashtekar and Varadarajan in the 2-1-1 dimensional setting We use here the rescaled scalar 
(f>s introduced above with the same dimensions of (j)g . As we will see they play symmetric roles in the final formulation 
of the model. The four-dimensional manifold where the previous action is defined has the form of a product Ai^ x Z 
where is a three-dimensional manifold orthogonal to the translational Killing vector = {d/dz)"- and Z = [zi, 22] 
is a closed interval in this direction (axis of symmetry). We have introduced also a fiducial metric g'^^^^ that provides 
us with an origin for the energy, ensures that the action is finite and fixes the asymptotic behavior of the fields in such 
a way that the Minkowski metric has zero energy. Finally h!-^l and h''^l° are the induced metrics on the boundary. 
Owing to the translation symmetry it is possible to rewrite the previous action as an equivalent one in 2-f 1 dimensions 
that can be interpreted as the Einstein-Hilbert action with two massless scalars after the conformal transformation 
gab = ^'^^Qab '^^ performed. Hence, our starting point to get the Hamiltonian will be 

^''' = T6^L, - \9''\'^<t>M<^<l^9)b~\9''Hd<l^s)a{d^^^^ {\h\^'^K-\h^\^'-'K^). (5) 

Here all the geometrical objects refer to the metric gab- The coupling constant G3 is the gravitational constant per 
unit length along the symmetry axis and in the following we choose units such that c = 1. In this three dimensional 
expression of the action, we notice that the scalar field term plays exactly the same role as the gravitational scalar. It 
is important to point out that both fields are coupled through the metric, but not directly (there are no cross terms 
in the action). 

To obtain the Hamiltonian we follow the procedure developed in 6] for the vacuum case. First of all, we choose a 
foliation of with timelike unit normal n", a radial unit vector f"^, and denote as a"' the azimutal, hypersurface 
orthogonal. Killing vector field (notice that this is not a unit vector). It is possible now to write the metric as 
gab = —na^b + rafb + -^CaCb (with R^ — gab<j"^cr^)- We also introduce two additional vector fields t°- and r°- defined 
as t°- = Nn"-+N'^f°' and = e^/^f°, where N is the lapse function, the radial shift, and at this point 7 is just an 
extra field (that will eventually coincide with the one introduced in the previous section). We impose the condition 
that the commutators of these new vector fields are zero, so we can define coordinates t, r, 9 and get the following 
consistency conditions 

d^N = d^N- - a.7 - ; [a, f]- = [a, n]"^ = ; 
n^drN + f'^idrN'' - dte'''^) + iVe''/^[f, n]" ^ 0. 

Finally the metric takes the form: 



gab = (iV^ - N^){dt)a{dt)b + 2e^/2^'-(di)(„(dr)b) -I- e^{dr)a{dr)b + R\da)a{da)b. 
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Now, if we take the boundary d^A^ to be orthogonal to the vectors f° and n°, it is straightforward to give an 
expression for the action in three dimensions in terms of the fields N , iV, 7, R, and (j)gy- 



8Lr3 Jti Jo V ^ 



R 
2N 
1 



27V 



dr dt 



AG 



t2 



3 Jti 



{Ne-'^/^R' -l)dt, 



where we have denoted dt with a dot and dr with a prime. We get now the Hamiltonian for the special case in which 
the boundary is taken to infinity (f — > 00) taking into account that the metric gat reduces to Minkowskian metric 
when iV = 1, iV'' = 0, 7 = 0, and R = assuming regularity in the axis, and the 2+1 dimensional asymptotic flatness 
conditions for the fields introduced in 0,113 ■ The Hamiltonian is then 



H 



7V'-e-7/2 



PrR' - 2p^ + P^fj' + (f)' Pg + <i)'^Ps 



r 1 



— (2i?" ~ 7'i?') - SGaPflp^ + -^p\ 

8Gr3 



^ ,12 , 2 



R 



R 



16G3 



dr 



^ 4G3^ ' 



where pn, p-y, pg, and ps are the momenta canonically conjugate to R, 7, <j)g, and (pg respectively and 700 
limf^oo It is easy to read both the constraints and the reduced Hamiltonian from the last expression. In order 
to proceed further we fix the gauge with the same conditions as in the absence of matter 

R(r) ~ r and P-y{r) = ; 

it is straightforward to show that these gauge fixing conditions are admissible. We can now solve the constraints to 
get 



7(7?) 



,2 ^ jSG^Pg)' ^ ^,2 ^ (8G3P.)- 



i?2 



i?2 



RdR, 



PR = ~Ps(l)'s -Pg^l^'g- 

Finally, the three-dimensional line element can be written as 



^dt^ + dR^] +R^da^ , 



(6) 



the reduced phase space is coordinatized by 4's{R), Ps{R), (t>g{R), and Pg{R); and the reduced Hamiltonian is 

1 



4G3^ 



where 



7o 



-7oo/2^ 

.2 , (8G3Pg)2 , .2 , (8G3Ps)2 



i?2 



i?2 



RdR. 



As we can see 700 is the Hamiltonian for a system of two free, axially symmetric scalar fields in 2+1 dimensions and 
the true Hamiltonian is a non-linear and bounded junction of this free Hamiltonian, similar to the one appearing 
in the absence of matter 12]. The Hamilton equations are 



-/2 Ps,g 

R 



Ps.g 



^e~^-/\R<^')' 



Though they are non-linear integro-differential equations they can be easily solved by realizing that 700 is a constant 
of motion. Taking this fact into account, we can perform a change in the time coordinate T = e.~^^l'^ t and rewrite 
them as 



hdr + dl + ^dn] 0.,^ = 0; 



(7) 
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describing two massless, axially symmetric scalar fields in 2+1 dimensional Minkowskian background. As we see this 
change in the time variable provides a one to one map from the solutions to a simple linear system to those of the 
non-linear one. This mapping encodes the non-trivial interaction present in the model. If we compare the classical 
evolution of two different sets of initial data we can see that they both correspond to the evolution defined by the 
free Hamiltonian 700 with times elapsing at different rates (defined by the conserved values of 7oo)- Notice that once 
we have a particular solution for equations (0 we have the freedom to "change coordinates" in the metric and write 
it in terms of t or in terms of T. This is not the case quantum mechanically, because the evolution of arbitrary 
states involves, in general, the superposition of Hilbert space vectors with energy (and time) dependent phases. This 
implies that quantum dynamics will be much more complicated than the classical one. In conclusion, the fact that 
the Hamiltonian is a function of a certain free Hamiltonian makes it both non trivial (we are, indeed, dealing with a 
coupled system) and solvable. 



IV. CANONICAL QUANTIZATION 

Once we have characterized the reduced phase space and obtained the classical Hamiltonian, we proceed to quantize 
the model. To this end we will use the Fock Hilbert spaces Tg^s associated with two different free, massless, axially 
symmetric scalar fields propagating in a Minkowskian background. These spaces are endowed with the usual creation 
and annihilation operators ag.s(fc), ^(fc) satisfying 

[ag{k), al{q)] = 6{k, q), [as{k), al{q)] = S{k, q); 

their corresponding vacua are denoted as |0)^''*. The Hilbert space of the interacting model is taken as the tensor 
product Ti — J-g ® Ts oi the Fock spaces corresponding to both scalars. We define the annihilation operators^ for 
modes of "gravitational" or "matter" types as Ag{k) :— ag{k) Is, As{k) :— Ig (g) as(fc), and the distribution-valued, 
field and momentum operators 0g.s(-R), Pg^s{R) 

poo 

^gAR) = V^G^ / MRk) [AgAk) + Al ,(fc)] dfc , 
Jo 

PaAR) = Y\f^,r ^^"^^^^ ^^'•'^''^ " ^'''^^^^ ' 

satisfying the usual commutation relations [(j)g^siRi),Pg.siR2)] = ih5{Ri, R2). Notice that we can construct states 
with a fixed number of quanta of "gravitational" or "scalar" type by acting with the corresponding creation operators 
on the vacuum state jil) = |0)f lO)** G 7i that is the minimum energy eigenstate of the quantum Hamiltonian'* 



H- 



4G, 



1-exp ( -AGzh j k [Al{k)Ag{k)+Al{k)As[k)] dfc) 



(8) 



This quantum Hamiltonian is a nonlinear and bounded function of the sum of the free Hamiltonians 

o 

fcit,,(fc)ig,,(fc)dfc 



-"0 



for two massless, cylindrically symmetric scalar fields in 2-\-l dimensions evolving in a fictitious Minkowskian back- 
ground. Their sum Hq + Hq is an observable but it is not the generator of the time evolution of the system. The 
physical evolution, from to to t, is generated by H and is given by the unitary evolution operator 



[/(t,to) = exp - 



i{t-to) 



1 



(9) 



This operator defines the S matrix of the system when we take the appropriate time limits. Its matrix elements on 
?i-particle states are straightforward to compute because these are eigenstates of the free Hamiltonian Hq = Hq + Hq. 



^ Creation operators are defined in an analogous way. 

* We have normal ordered the exponent because, otherwise, the Hamiltonian is trivial. 
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As we can see the only matrix elements -involving state vectors with a definite number of both type of quanta- that 
are non zero are those connecting states with the same number of particles of each type; hence there is no conversion of 
quanta of one type into the other. At this point, it is important to reflect upon the interpretation of these elementary 
excitations of the fields. One must be careful, for example, when interpreting states such as |0)^ (8) |$)* because one 
should not be lead to think of them as matter (represented by |$)'*) evolving in a certain background (given by |0)^). 
The system that we are considering is a coupled one, and hence one is not entitled to think of the metric and matter 
fields as independent objects. They are coupled in the equations and, in particular, the classical metric depends on 
both the gravitational and matter fields. Conversely, the evolution of the scalar field depends on the metric. If we 
want to approximate the Minkowski metric, we have that the state that most closely resembles it is the vacuum of 
our total Hilbert space jfJ) [by the way, this is the only coherent state of the system that we know under the evolution 
©]. In the next sections we will make extensive use of states of the type |0)^ (8) \k)'^ consisting of tensor products of 
the vacuum state of one of the Fock spaces and a one particle state of the other. As we will show in the last sections 
of the paper these single particle states are the closest ones to Minkowski if one wants to incorporate an extra element 
that can be used to explore and describe the geometry of the quantized model in an operational way. 

In the following it will be convenient to explicitly keep the length scale of the system G — G^h in the mathematical 
expressions of the relevant objects. We will nevertheless use units such as fi, = 1. With the time evolution (|5J defined 
by the Hamiltonian (jHl the annihilation and creation operators in the Heisenberg picture are 



As,g{k), 



it (k- 1, to) = UHt, to)As,g{k)U{t, to) = Al (k) exp i{t - to)E{k)e~^''"° 



where 



Eik) :=^(l-e- 



~4G/c 



and Hq := 



• 



Then the scalar field operators that describe the gravitational and the massless scalar field degrees of freedom at time 
t are 



s,g{R;t,to) = V4G / Jo(i?fc) 



As,g{k;t,to) 



Mjk;t,to) 



dk. 



V. TWO-POINT FUNCTIONS 

As commented in the introduction one of the main goals of the paper is to find a way to recover a physical picture 
of spacetime in the quantized symmetry reduction of gravity coupled to matter discussed above. We want to find 
ways to describe a quantized spacetime geometry in an operational way much in the same way as one explores a 
classical spacetime geometry by using test particles. To this end it is useful to have objects playing the role of 
particle propagators or, even better, one-particle states that could allow us to define suitable wave functions with a 
straightforward interpretation as spatial probability amplitudes. Hopefully the time evolution of these objects may 
give us some idea about the physical effects of quantizing the gravitational field and also tell us something about 
how the classical macroscopic geometry emerges. As the reader may expect this is not easy; in fact the problem of 
finding suitable position eigenstates in the usual Minkowskian QFT's is already non-trivial. Our strategy will be to 
use some of the objects introduced in the discussion of these issues in the traditional approaches to QFT (such as 
two-point functions or Newton- Wigner states), interpret them in our framework, and use them to obtain a physical 
picture of the quantized geometry and gravity. Also, in order to disentangle genuine quantum gravitational phenomena 
from artifacts introduced by the symmetry of the problem and the reduction process we will compare the relevant 
objects to the corresponding ones in a model of a quantized, axially symmetric, massless free scalar field moving in 
a Minkowskian background. This comparison is conceptually simpler if one works both with gravity and the scalar 
field. 



The vacuum expectation values of the product of two fields at different spacetime points can be interpreted, at 
least in an approximate sense, as propagation amplitudes for particles or field quanta created at a certain event to be 
found at another. This is so because the (Schrodinger picture) scalar field operators that describe the gravitational 
and scalar degrees of freedom are 
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(i?) V4G / Jo (Rk) \As,g (k) + (fc) 



dfc 



and their action on the vacuum \n) = |0)^ |0)'* satisfies 

1 



/AG 



<^,,g(i?)|fl) = / dk JoiRk)Al Jk)\n) ^ / dkJoiRk)\k) 



^0 



These are hnear superpositions of (ortlionormal) states \k)s^g ■— A'g ^{k)\n) with well-defined "radial momentum" k. 
Notice that 

MRk) = ,Jk\$s,g{Rm) (10) 



is a solution of the radial part of the Schrodinger equations for states with zero angular momentum in two dimensions 

[dl + ^dn + k^]MRk) = 0. 

K 

If we consider a small volume element /S.V at a distance R from the symmetry axis the value of J^{Rk)AV is 
proportional to the probability of finding a particle of type s or g inside it. Notice that this is not in general the 
probability to find the particle in a thin cylindrical shell of radius R. 

In order to consider the amplitude for propagation of quanta of the matter scalar we could consider more general 
situations, for example we might take states of the form |c)^(8) with |c)^ a suitable "coherent" gravitational state 
under the evolution defined by the dynamics of the system. Notice, however, that the physical interpretation of such a 
state is not completely clear and, in particular, one should not be led to think that the gravitational part of the state 
(say |c)^) fixes the geometry (or a suitable classical approximation thereof) and the matter part the scalar field; in fact 
both parts of the state vector contribute to fix the metric on one hand and the matter field on the other (classically 
this can be understood by realizing that the metric depends on both the gravitational and matter scalars). 

In the following we will consider the case ^2 > ti and interpret the matrix element (r2|(/)s,g(i?2; ^2, ^o)'/'s.g(^i; ii, ^o)|^^) 
as the (approximate) probability amplitude of a particle created at a point at a distance Ri from the axis in the 
instant of time ti to be detected at a another point at i?2 distance in the instant of time ^2- We can now obtain in a 
straightforward way^ 



{n\^s,g{R2;t2,to)A,g{Ri;ti,to)\n) = AG / JoiRik)Jo{R2k){n\A,,g{k;t2,t„)Alg{k;ti,to)\n)dk 

Jo 

/•oo 

= 4G / Jo{Rik)Jo{R2k)exp[~i{t2-ti)E{k)]dk. (11) 
Jo 

We have introduced an initial time to that will not appear in the final expressions of the matrix elements that we will 
consider here so that in the following we will write {i}\(j)s,giR2,t2)4>s,g{Ri,ti)\il) . 

Let us consider the integral (|11() . First of all it must be said that it is not possible to compute it in closed form 
although there are suitable ways to compute it numerically and approximate it by means of asymptotic expansions 
jsj. The relevant parameters in the integral are: i?2, ^2 — ti -the arguments of the two-point function-, and AG 
that sets the length scale.® In view of this it appears to be appropriate to refer both length and time to this scale and 
introduce the adimensional variables Pi = ^ , P2 = ^ , and t = together with the change of variables q = AGk 

that gives a dimensionless integration variable. In this way we can rewrite Hll|l as 

/>CXD 

{n\k,g{R2,t2)$s.g{Ri,ti)\n) = / Jo(pig)Jo(p2g)exp[-iT(l-e-'')]dg, (12) 

Jo 



The creation and annihilation operators as written above are specially suitable for this computation and similar ones where the states 
are eigenstates of the Hamiltonian operator. 

In fact 4G plays the role of the Planck length. Notice also that it also sets the time and energy scales as we are taking units such that 
h=c=l. 
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where pi and p2 are to be considered as functions of i?i and R2 as defined above. An interesting consequence of H12() 
is that we can obtain the vacuum expectation value of the commutator of Heisenberg picture field operators by taking 
its imaginary part (as we did in to discuss the microcausality in this system). 

One can consider, in principle, the numerical computation of this type of improper integral but this is rather difficult 
due to the oscillating nature of the integrand. In spite of this there arc efficient ways to do it, as the one described 
in the appendix, by rewriting it as an integral over a torus plus a rapidly convergent improper integral. On the other 
hand, the advantage of using asymptotic approximations in some relevant parameters lies on the fact that it gives 
the limiting behavior at large scales in some physically relevant regimes and also allows to get numerical estimates 
in a rapid manner. It is possible to consider separate expansions for each of the parameters pi, p2, and r which are 
valuable in the sense that it is possible to study the behavior of the two-point function when only one of them is 
taken to be large. It is also possible to write down an expansion in which all of them are simultaneously large while 
keeping their relative values. Here this approximation corresponds to large length and time intervals as compared to 
the scale set by G. Let us consider these separately: 



A. Asymptotic expansions in pi or p2 



For large values of pi and p2 the asymptotic behaviors of 1)12(1 are respectively 



1 

Pi 



pI 



IT 

y 



+ 0{pl^) and 



1 

P2 



IT 



Pi 



'p1+t--- 
4 2 



+ 0(p2-6), 



they can be obtained in a straightforward way by considering H12|) as a standard /i-transform with asymptotic param- 
eters pi or p2 and using Mellin transform techniques |l3j | . The imaginary part of the previous expressions corresponds 
to the vacuum expectation value of the field commutator. It does not show the sharp discontinuity present when one 
considers axially symmetric massless scalar fields in a Minkowskian background and, hence, is a quantitative measure 
of the spreading of the light cones expected in a quantized theory of gravity as discussed at length in |^ . 



B. Asymptotic expansions in r 



In order to study the asymptotic behavior in r for \12\ it is convenient to consider two separate cases: Both pi and 
P2 different from zero or only one of them equal to zero (the other must be different from zero because, otherwise, the 
integral is badly divergent). In the first case the expansion can be obtained by first writing H12II as an /i-transform, use 
the Mellin-Parseval formula to obtain a representation for it as a complex contour integral, and split it (by displacing 
the integration path close to the branch cuts of the integrand) into several integrals that can be studied by standard 
Mellin-transform techniques. The details of this procedure -that allows to obtain the asymptotic behavior to any 
order- can be found in (3j. The result in the present case is 



1 



27ry'pip2 logr 



(exp (|(P1 +P2) - + T- (Pl +P2)l0gT])r[-i(pi +P2)] 

/ TT TT \ 

+ exp(^--(pi +P2) -r- (pi +p2)logr] jT[i{pi+ P2)] 

+ exp(^(pi -P2) -«[t+ (p2 -pi)logT]^r[i(p2 -pi)] 
+ exp (|(p2 -pi)-i[T+ (pi - P2) logT])r[z(pi - P2)]} + 0(1/ log' t). 



When either pi or p2 are equal to zero we respectively obtain 
1 



A/27rp2 log T 
and 

1 

V27rpi log T 



exp |^-P2 - + r- p2 logT)Jr(-ip2) + cxp |^--p2 + «(- - r - p2 logT)J r(ip2)| + 0(1/ log^/' t) 



exp 



:Pi 



T- Pl log r) 



r(-zpi) +exp [-|pi -r-pi logr)] r(zpi)} -1-0(1/ log^/'r). 



As we can see for fixed values of pi and p2 the decay when one of them is zero is slower (an inverse power of v'log r) 
than the decay when both pi and p2 are different from zero (an inverse power of log r) . We can compare this result 
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with the one corresponding to a massless, axially symmetric free scalar field'^ 

/"OO 

(O|0(i?2,i2)'^(i?l,tl)|O) = / Mpiq)Mp2q)eM-iTq)Aq. 

Jo 

For values of r satisfying t > pi + p2 this integral is equal to^ 



2« r^/ / 4pip2 




TTv^T-^ - (pi - Pif vV ^ ~ 



and for large values of r it behaves asymptotically as — i/r, i.e. it falls off to zero much faster than H12fl . As a 
consequence of this we interpret the very slow decay of the two-point function for quantized Einstein-Rosen waves 
as an enhanced probability amplitude to find quanta (either of the gravitational or matter field) in the symmetry 
axis. This is a gravitational effect as it is no present for a quantized, axially symmetric, massless scalar field in a 
Minkowskian background. 



C. Asymptotic expansions for pi, p2, and r simultaneously large 



We discuss now the obtention of an asymptotic approximation that is valid in a "gravitational classical limit" 
corresponding to taking pi, p2, and r large while keeping their relative values. This is equivalent to consider values 
for i?i , i?2 , and t2 — 1\ which are much larger than the "Planck scale" provided by AG; it is in this sense that we talk 
of a classical limit here. To this end let us rewrite the integral in H12I) as 



Jo(pig)Jo(p2g)exp[-iT(l - e '^)] dq 



Jo{Xriq)jQ{Xr2q) exp[-«Ai(l - e '')] dq 



where pi = Xri, p2 — Ar2, and t ~ Xt with r2, t fixed, and A taken as a new parameter that we will consider large 
(we will use it as an asymptotic parameter A 00). The last integral can be written as 



47r2 



dg i) dzi i) dz2 
J 72 



cxp A —(21 ) + -Tr{z2 )+ite « 

Z1Z2 V L 2 zi 2 Z2 



f dz 


■X , 






— z - 
L2^ 





by using the well-know representation of the Bessel functions as contour integrals (over contours 7 that enclose z — 0) 

As discussed in it is useful to chose the contours 71.2 in the complex plane region satisfying 3fi(z — 1/z) < 0, z € C. 
A suitable asymptotic expansion for it in A can then be obtained by following the procedure outlined in We 
only quote the result here after reabsorbing the parameter A (and, hence, expressing the integral again in terms of 
pi, p2, and r). To do this we need to define three different regions that cover the (pi,p2,r) space, ^ referred to in 
the following as I, II, and III, and defined by the conditions: r < |p2 — pi| for region I, |p2 — pi| < r < pi -(- p2 for 
region II, and t > pi -I- p2 for region III. The asymptotic expansion in each of them is given by the sum of a leading 
contribution (~ 1/A) and a first order correction behaving asymptotically as 1/A'^/^ or 1/A^. The leading contribution 



Here we denote the vacuum state of the system as |0). 

/•7r/2 

Complete elliptic integrals of the first, second, and third kind are respectively defined as K{k) = / — , E(k) = 

Jo V 1 - fc2 sin^ e 

/'7r/2 /'7r/2 

/ VI - fc2sin2 6»de, and n{n|fc) = / - . 

-'0 Jo {1-n sin^ 0) V f - fc2 sin^ e 

We restrict ourselves to positive values of t -the extension to negative values is straightforward- and avoid the boundaries between 
these regions. 
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in the different regions is 



Region I : 
Region II : 

Region III : 

and the first asymptotic correction 

'2v/(pi+P2F 




(13a) 
(13b) 
(13c) 



Region I : 



2piT' - 3r4 + 2piT 



,2^2 



(Pl + P2 - T)^iPl - P2+ t)^(-Pi + P2+ tY{PI + P2+ t)^ 

2r2 



ViPi+P2y-r^pt + ir' Ply - '^pUpI + r^)] 




Region II 



Region III 



g4f-r + |p.-pi|(l+log^;jf^)] 
^/2'KpiP2\pi - P2\ log 
1 



\P2-Pl\ 
T + |P2-pl|(l+log 



^)] -4f+r-(pi+p.)(l+log^)] 



V27rpip2 



y/\pi - P2|log 



|P2-P1 



+ /92 log 



(14a) 
(14b) 
(14c) 



P1+P2 



The imaginarypart of these expressions gives the vacuum to vacuum matrix elements of the field commutator already 
discussed in Also it is important to point out that (|13|l is, precisely, the two-point function for an axially 
symmetric massless scalar field evolving in a 2+1 dimensional Minkowskian background. In fact, this is the leading 
contribution to the two-point function -corresponding to an ordinary quantum field theory for a massless scalar field 
in a Minkowskian background- whereas (|14|) provide the first asymptotic corrections. 

It is also possible to obtain similar expansions when either pi or p2 are zero. In this case we only have a single Jq 
function in the integrand and the computations are greatly simplified. The result is 



e{p-T) 



lT{p^ + 2t2) 



Vp^ 



(p2 



2)5/2 



+ Oir - p) 



exp(i[plog^ -T + p]) 



log^ 



(15) 



where p is the remaining non-zero radial parameter. Notice that this cannot be obtained by simply putting pi = or 
P2 = in CTTl) . 

Figures ^ and [21 show the behavior of the two-point function. One can compare the exact values obtained by 
numerical computation and the approximation given by the asymptotic expansions. We also compare its value with 
the one corresponding to an axially symmetric massless scalar field evolving in a Minkowskian background. We 
emphasize again here that the value of these asymptotic approximations relies on the fact that they give the exact 
behavior of the two-point functions in the relevant limits. The most important physical information that can be 
gleaned from these plots is a significant enhancement of the probability to find field quanta (cither gravitational or 
matter) in the vicinity of the symmetry axis (defined by p = 0) as compared with the result for an axially symmetric 
massless scalar field in a Minkowskian background. This is specially remarkable because far from the axis (both pi and 
P2 large) the dominant contribution to the two-point function is given by the one corresponding to the free massless 
field. The asymptotic analysis for p = Q shows that the symmetry axis is exceptional in the sense that one does not 
recover the free- field result in the asymptotic limit. This was already noticed in previous studies of the microcausality 
of the system; in this context the interpretation of this phenomenon is the enhanced probability mentioned above. 
Another feature that stands out in the figures is the singularity at pi = P2- This is expected on general grounds as 
a generic behavior in quantum field theory. It can be removed by introducing suitable regulators (see f^). Here we 
can, in practice, identify and isolate the distributional behavior of the relevant objects at these points so we will not 
introduce regulators explicitly. 



Other interesting two-point functions that we may study are mixed ones involving both the gravitational and scalar 
fields. We then consider matrix elements of the form {f\(t)s{R2\ t2,to)(t)g{Ri;ti, to)\n) for some state |/) different from 
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FIG. 1: This figure shows the absolute value squared for the two-point function for fixed values of p2 and r in terms of pi. It 
shows the approximate probability of finding a field quantum at a small volume centered around pi after a certain time lapse r 
if its position was p2 at r = 0. Here radial distances and time are measured in units of 4G. We also compare the exact values 
obtained by numerical computation (labelled as "Num."), the approximation given by the asymptotic expansions (labelled as 
"Asymp."), and the one corresponding to an axially symmetric massless scalar field evolving in a Minkowskian background 
(labelled as "Free"). The most salient feature is the significant enhancement of the probability for pi = 0. The dot on the 
vertical axis corresponds to the value of the two-point function on the axis pi = 0. The remarkable quality of the approximation 
provided by the asymptotic expansions given in the paper can also be seen except in the boundary between regions in the pi, 
p2, and r space, where the asymptotic expansions are divergent as expected on general grounds. The inset shows in detail the 
comparison of the exact values of the two-point function and the asymptotic approximation given by I13I14II . It is worthwhile 
to point out that even though one does not expect the approximation to be valid for small values of pi or p2 its general behavior 
is well described by it. Notice also the expected singularity at pi — p2- 



the vacuum. The simplest non-zero matrix elements are obtained by choosing 

/•oo poo 

1/)=/ dkj dkgf{k,,kg)Al{ks)Al{kg)\n) 

Jo Jo 

satisfying the normalization condition 

/•oo poo 

/ / \f{ks,kg)\^dkgdks = l. 

Jo Jo 



Diagonal matrix elements of this type are always zero -for both scalar and gravitational modes— because they involve products of an 
even number of creation operators with an odd number of annihilation operators (or viceversa) that act on the vacuum to give zero. 
Also expectation values between the vacuum and one particle states of either type are easily seen to vanish. 
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25 50 75 100 125 150 



FIG. 2: This figure shows the absolute value squared for the two-point function for a fixed value of p2 and pi = 0. The most 
interesting feature is the fact that the squared amplitude does not behave in mean as the one corresponding to the "free part" 
as was apparent in the previous picture. The dot at r = 150 corresponds to the value of the two-point function at the symmetry 
axis shown in figure Q Notice the different scale used in both plots. 



They can be computed easily to give 



{f\MR2;h,to)$g{Ri;h,to)\n) 



4^ J^'^'^''J^'^'^3MP2qs)Mpiqg)f{-^,-^)exp[i[T2o{l~ e «-')e + tio{1 - e 



(16) 



where we have mtroduced adimensional variables and parameters as above: kg = j^, kg — j^, pi — p2 = j^, 



no 



= *° , and T20 = ■ It is also interesting to write the matrix element of the commutator 



{f\[MR2;t2,to),MRi->ti,tom) 

= 4G / dfc, / dkg J0{R2k,)j0{Rlkg)f{ks, kg) 



1 



"'0 

OO pOQ 



gi[(t2-to)B(fc.)e-'"='=9+(ti-to)B(fcg)] _ gi[(t2-to)B(fe.) + (ti-to)S(feg)e~'"^'=''] 



Jo 



qs qg 



= ^ lMsjdqgjQ{p2qs)Mpiqg)f{^,-^) 



4G' AG' 



gi[r2o(l-e-''=)e-99+rio(l-e-99)] gi [r2o(l-e-«= )-|-tio (l-fi-'s )e-9= ] 



(17) 



Some conclusions can be reached by considering some simple choices for the function / although a more detailed 
analysis would require us to get suitable asymptotic expansions in terms of a general /. To this end let us consider 
the normalized function 



f{ks,kg) ~ -^X[ko^-k/2,ko,+k/2](^s)X[kgg-k/2,kog + k/2](^g) 
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where xv is the characteristic function of the set V and fc is a constant with dimensions of inverse length. For values 
of k small enough that 

MR2ks)MRikg) exp (^i[{t2 - to)E{ks)e-*^''^ + {h - tQ)E{kg)] 
is essentially constant in the effective integration region the value of the matrix elements (|16I17|) are respectively 

^GkMR2kos)MRikog){e'^^''^'''^^^''^''^'''^^^^ 
and their squared amplitudes 

8G2P j2(i?2A:o.) Jo'(i?iA:og){l - cos[4G(t2 " h)E[kas)E{k,,g)]}. 

The first important point to notice here is the fact that these quantities, and in particular the commutator, are 
generically different from zero. This is an additional indication of the fact that we are dealing with an interacting 
(i.e. non-free) theory. It is also easy to see that the commutator is zero when ti = t2 because 

exp (i{t2 - to)[E{ks)e-~'^^^o + E{kg)]] = exp ( ^(^2 - to)E{ks + kg)] = exp ( ^(^2 - to)[E{kg)e-'^^^' + E{k. 



Finally it is remarkable that these functions do not display the type of causal behavior that we have found for the 
commutator of fields of the same type (by looking, for example, at the vacuum to vacuum expectation value). Their 
magnitude is of order G^k^ so, at least in this approximation, are small in the scale set by the natural length scales 
of the model. 



VI. NEWTON-WIGNER STATES: PROPAGATORS 



The main drawback of the two-point functions discussed in detail in the previous section is the fact that their 
interpretation as probability amplitudes is only approximate. The reason behind this is the fact that the state vectors 

0,,g(i?;t,io)|f^) 

do not constitute an orthonormal set. This is at the root of the well known problem of localization in relativistic 
quantum field theory that has been solved in a more or less satisfactory way by the so called Newton- Wigner states 
[3 • These constitute an orthonormal basis of position eigenstates for a certain choice of an inertial reference system 
in a Minkowskian spacetime. The most important interpretive difficulty with them is the fact that they cease to be 
localized under Lorentz boosts Q but the assumptions upon which their construction is based are so natural that it 
is difficult to believe that a better solution to this problem might exist. 

The purpose of this section is to build localized states analogous to the Newton- Wigner ones for our model, that we 
will label here as | i?) . We have now a reduced spacetime symmetry group so the problem of the behavior of localized 
states under spacetime symmetry transformations is partially alleviated. The physical interpretation of propagation 
amplitudes built with this type of states is clear in the sense that they are now proper probability amplitudes. Our 
point of view here is that the discussion of the two-point function (f2|(/)s(i?2, i2)0s(-Ri, ii)|f^) together with the matrix 
elements (i?2|f^(i2, ii)|-Ri) can give us relevant and robust information about the motion of field quanta in position 
space and provide meaningful information about the transition between quantum and classical geometry in quantum 
gravity. It should be pointed out here that the availability of a position space orthonormal basis allows us to define 
position state normalized wave functions 

1*00 /"OO 

1*)=/ <iR'^[R)\R), = «'(i?), / |*(i?)pdi?=l 

Ja Jo 

and study their time evolution (in the Schrodinger picture) given by |5'(i)) = U{t,to)\'i>{to)). 
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Here we will follow a simple procedure inspired in . Let us write 

/■oo 

\R) = / dkfik)MkR)\k) 



Jo 

where \k) := |0)^ (8) are scalar matter field quanta and we have made use of the fact that Jo{kR) is a solution of 
the radial part of the Schrodinger equation for states with zero angular momentum in 2 dimensions: 

[dji + ^dR + k^]Jo{kR) = 0. 

Once we make this choice the function f(k) is fixed by the orthogonality condition (i?2|-Ri) = '5(i?2, Ri) which implies 
|/(fc)p = kR, and hence f{k) = \/kRe^'^^''\ Without loss of generality we will take iy{k) = 0. We finally get then 



= / dk{kR)^^^Jo{kR)\k). 
Jo 

In the following we will study matrix elements of the form 

/>oo 

{R2\U{t2,ti)\Ri) = ^/RiR2 / kJo{kRi)Jo{kR2)exp[-i{t2-ti)E{k)]dk 

Jo 

= J%Mpiq)Mp2q)eM-ir{l-e-'^)]dq (18) 



^^S{pi,P2) + '^Q^ ^ " qJo{piq)Jo{p2q) exp(ire - 1 dq 



The last integral in the previous expression converges very quickly as the integrand has an exponential decay. Notice 
also the singularity at pi = p2 that was also present for the two-point functions. As happened before we cannot give 
a closed form expression for the integral (|18|l although, again, it can be computed numerically by essentially the same 
methods used in the previous section. We can also obtain asymptotic approximations of the types discussed above 
that help us understand precisely the behavior in several important physical regimes and in the "classical limit" . 

An important question is the meaning of this probability amplitude. In the case of the two-point function we 
highlighted the interpretation of (|10() as the radial part of a wave function with zero angular momentum for a free 
two-dimensional particle. Now^^ |(fc|i?)p = kRJg{kR) and the appearance of the R factor suggests that the correct 
interpretation for the amplitudes given Newton- Wigner states is that they describe the probability to find field quanta 
inside thin cylindrical shells at a distance R. This means that we will have to introduce appropriate factors of Ri and 
i?2 to compare the two-point functions of the previous section with {R2\U{t2,ti)\Ri) . Specifically we will study 



^=^{R2\U{t2M)\Ri) = / <zJo(pi<z)Jo(p2g)exp[-ir(l-e-?)]d(?. 
V-Ki-n2 Jo 



(19) 



In the following we give the different asymptotic expansions for this propagator in the same regimes described in the 
previous section. 



A. Asymptotic expansions in pi or p2 

For large values of pi and p2 the asymptotic behaviors of H19|) are respectively 

obtained, again, by a straightforward application of Mellin transform techniques. When this is compared to the 
asymptotic behavior obtained for the two-point function we see that they qualitatively agree for the imaginary part 
but they are quite different for the absolute value or the real part. This is not unexpected as the interpretation of the 
two-point function as a probability amplitude is approximate. 



i + 



Pt 



I T 

'6 ^ 1~ ^ 2 



IT 



p\ 



Notice that (fe|-R) is not a solution to the Shrodinger equation. 
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B. Asymptotic expansions in r 

As before it is convenient to study separately the case in which both pi and p2 are different from zero and the one 
in which either pi or p2 are zero. By using the same methods as before we find 



{exp + P2) - + T - (pi + P2) logT]^r[-i{pi + P2)] 

~-^{Pi + P2) + - T - (pi +p2)logT] jr[i{pi + P2)] 
+ exp(|(pi -P2) (P2 -pi)logT]^r[i(p2 -Pi)] 

+ exp(|(p2 - pi)-i[T+{pi - p2)logT]jT[i{pi -P2)]} +0(1/ log r). 

When either pi or p2 are equal to zero we respectively obtain 



/logr / 

V ^ 



rP2 



( J + r - p2 l0gT)]r(-ip2) + exp [-^P2 +l(^-T-p2 logT)]r(ip2)} + 0(1/ log^/' t) 



and 



' logr 
27rpi 



exp 



j(- + T - pi logr) r(-ipi)+exp 



2^1+^(4 



T - Pi log r) 



X^Pl)} + 0(l/log 



1/2, 



As we can see for fixed values of pi and p2 when one of them is zero the asymptotic behavior of (|19l) for large values 
of T consists of an oscillating function times VIogT, also it must be noted that the oscillating part is precisely the one 
that appeared in the r-asymptotics of the two-point function. For pi and p2 both different from zero the asymptotic 
behavior is given by factors that are purely oscillatory in r. Again these coincide with those that appeared in the 
study of the two-point function. We see now that the value in the axis grows (very slowly) and it has a constant 
amplitude everywhere else. We see again that there is an enhancement of the probability to find the particle in the 
axis relative to the rest of the values of pi and p2 . We can compare this result with the one corresponding to the free 
axially symmetric scalar field in 2-1-1 dimensions given (for r > pi -I- P2) by 



2r 



^[(pi + P2)" - t2] - (pi - p2)^ 




4piP2 



(pi - P2Y 



As we can see for large values of r it falls off to zero as — 1/r^, much faster than the asymptotic expansions that we 
have already found. Again we interpret this result as an enhanced probability to find field quanta in the vicinity of 
the symmetry axis. We see that although the analytic expressions that we have obtained are different from the ones 
corresponding to the two-point function the qualitative conclusions regarding the behavior in the vicinity of the axis 
are the same. This strongly suggests that we are seeing a genuine quantum gravitational effect. The most significant 
difference between both types of results is the fact that the probability at the axis decays very slowly if one considers 
the two-point function whereas it slowly grows if one uses the Newton- Wigner propagator. We do not perceive a 
contradiction here because of the approximate interpretation of the result for the two-point function and the necessity 
to look at probabilities over spacetime regions. In fact we will look at this again in the next section when we consider 
the evolution of actual wave functions. 



C. Asymptotic expansions for pi, p2, and r simultaneously large 

We discuss now the obtention of an asymptotic approximation that is valid in a "classical limit" corresponding to 
taking pi, p2, and r large while keeping their relative values as we did for the two-point function. To this end let us 
rewrite the integral in (|19|l as 

/•oo 1*00 

I qjQ{piq)jQ{p2q)eii-p[~iT{l- e^'i)]dq= I qJa{\riq)jQ{\r2q) Gyip[-i)^t{l - e^'^)]dq 
Jo Jo 
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where pi = Ari , p2 = Ar2 , and t ^ Xt with ri , r2 , t fixed, and A taken as a new parameter that we will consider large 
(we will use it as an asymptotic parameter A — > oo). The last integral can be written as 



itX 



47r2 



dq (p dzi d) dz2 

Jfi J 12 



q ( I 1 qr2, 1 , 
exp ( A — (zi ) + -7r[Z2 ) 

ZiZ2 \ I 2 Z\ Z Z2 



ite 



by using, again, the representation of the Bessel functions as contour integrals. The leading contribution in the 
different regions is 



Region I : 
Region II 
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7r[(pi - P2 f - t2]v^(pi +P2) 
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4pip2 
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Region III 
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7r[(pi + P2)2 - r2] - (pi - p2? 
and the first asymptotic correction 
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Region I : 
Region II : 
Region III 
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^27rpip2|pi - P2I 

1 I ^ii-r+\P2-pi\{l+ioej^^)] ^-.,[f+r-(pi+p2)(l+log^)] 
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The previous expression for region I can be explicitly written in terms of complete elliptic integrals of the first and 
second kind with coefficients that are square roots of rational functions of pi, p2, and t; as they are rather lengthy 
we do not give them here. As before the leading contribution corresponds to the Newton- Wigner propagator for a 
massless axially symmetric scalar field evolving in a Minkowskian background. It is also possible to give asymptotic 
expansions in the case when either pi or p2 are zero. They are 



0{p-r) 



(p2 



3t^3p^ + 2t^) 

^)3/2 ^ 2(p2 -r2)7/2 



+ Sir - p) 



(^2 _ ^2)3/2 p 



^e'(''-^+''i°s^\ Aog^ 



where, as above, p is the remaining non-zero radial parameter. Notice that this last expansion cannot be obtained by 
simply putting pi = or p2 = in I|20I21|I . 

The most interesting feature of the propagator is its behavior in the axis and at pi — P2- Let us consider first the 
behavior at pi = P2 . Here we find the singularity expected on general grounds due to the orthonormality property of 
the Newton- Wigner vectors. This can be identified as the delta function appearing in (|18ll . In addition to this we see 
a clear tendency of the probability amplitude to remain concentrated around the region pi = p2. We interpret this as 
an effect of self gravity that tends to favor the concentration of matter. The spreading of the amplitude, on the other 
hand, can be interpreted as quantum mechanical diffusion similar, but less extreme, than the familiar one for particles 
in ordinary quantum mechanics. At the axis pi = (or p2 = 0) we see that once the amplitude grows as a consequence 
of the gravitational collapse of the initial matter distribution there is a tendency to have a large probability to find field 
quanta at that position (in fact the amplitude grows!). To really assess how the probability of finding particles near 
the symmetry axis evolves we will consider in the next section the evolution of proper normalizable wave functions and 
confirm, indeed that the probability is enhanced but decays very slowly in time. These effects bear some resemblance 
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FIG. 3: This figure shows the square modulus of the Newton- Wigner propagator (divided by factors of -Ri,2 introduced in 
order to compare it with the two-point function discussed in figure^!. We can see several interesting features: 1) An enhanced 
amplitude at the axis similar to the one already seen for the two-point function, 2) A large amplitude at pi = P2, even when 
the delta function at this position is subtracted; we interpret this a self gravity effect in a region of hight matter density, 3) 
Beyond the position corresponding to pi = 200 the amplitude decays very quickly; this marks the position of the light cone. 
Notice that even though the amplitudes in the "free case" -corresponding to the propagation of a massless axially symmetric 
scalar in a Minkowski background- diverge in some regions in the (p, r) plane they remain finite in our quantum gravity model. 



with what one expects to find in the study of gravitational collapse and black hole evaporation.^^ It is important 
to compare the results of this section with the ones derived in the study of the two-point function and the free 
axially symmetric massless scalar field. Regarding the first it should be pointed out that the qualitative agreement 
between the two pictures is very good. We see in both cases the enhancement of the probability at the axis, the 
singular structure in the vicinity of pi = P2, and the motion along null radial geodesies (manifest itself as a significant 
probability to find particles at the "classical" light cone). When the results are compared to the "free" massless case 
the gravitational phenomena that we have interpreted as an increased probability due to self-gravity are conspicuously 
absent but those related to the causal structure and microcausality of the system (i.e. enhanced probabilities on the 
light cones) are still present. 



The main difference is the apparent absence of a horizon or something behaving, at least in an approximate way, as one. 
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FIG. 4: This figure siiows tlie square modulus of tlie Newton- Wigner propagator (divided by factors of -Ri,2) at the axis. The 
dot corresponds to the value shown in figure 



VII. TIME EVOLUTION OF A (RADIAL) POSITION STATE WAVE FUNCTION 



For superpositions of one particle states the time evolution of the wave function = {R\'i'{t)) is given by 

/>oc />oo 

^{R,t) = {R\U{t,to)\^{to)) = / / k\/RRJo{kR)Jo{kR)exp[-i{t-to)E{k)]^{R,to)dRdk. 



'0 



In the following we want to study the evolution of one particle wave functions of this type. This will be, in fact, the 
quantum test particles that we will use to describe the quantum geometry of our spacetime model. In principle we 
can make any choice of initial wave function ^'(i?,to)- However, in order to obtain closed expressions for the wave 
function (or at least as simple as possible) we will concentrate on a specific choice that satisfies several reasonable 
requirements: The possibility of having some control on the position of the peak of the probability distribution, the 
possibility of controlling the width of the wave packet, and the possibility of performing (some) integrations to get a 
manageable closed form for it. A (normalized) function satisfying these conditions at to is 



/ 2R 

*(i?,to) = \ — 2 ^['■i.'-2](-^) (withr2 > ri) 



that gives 



/ 2R 

^{R,t)^J- 2 / Jo{kR)[r2Mkr2)-nJi{kri)]exp[-i{t~to)E{k)]dk. 

V ^2 ^1 Jo 

It is convenient to rewrite this in terms of adimensional objects by introducing the following redefinitions p = ^ 



<^2 = jh, en = jh, T = K^T^ and the change of integration variable q = AGk. Writing t) := ^{AGp, AGt + to) we 
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finally get 



4G( 



2p f°° 

— 2T / -'^o('7/o)[o-2Jl('7cr2) - CTi Ji(qcri)] exp[-iT(l - e '*)]dq 

^11 Jo 



(22) 



that satisfies, for all r, the normalization condition 



4G / \^{p,T)\'dp=l. 
Jo 

It is straightforward to see that when r = we recover the initial wave function at t = io- Our goal now is to 
extract information about the evolution of this wave function. In particular we are interested in the behavior near 
the symmetry axis to see if the enhanced probability suggested by the analysis of the two-point functions and the 
Newton- Wigner propagator is present. We also want to find out if the evolution of this wave function somehow defines 
a classical trajectory in the spacetime that can be used to define in an approximate way a physical notion of geodesic. 
Finally we want to compare the result with the one obtained for the axially symmetric massless scalar that we are 
using to disentangle quantum gravitational effects from more prosaic behaviors. As in the previous sections we will 
use asymptotic approximations in different regimes to extract the analytic behavior of ip{p-, r). 



A. Asymptotic expansion in p 

For large values of p the asymptotic behavior of H22|l is 

9 



4G p3/2 



2p2 



1 IT o , 



+ oip-'/') 



obtained, again, by a straightforward application of Mellin transform techniques. When this is compared to the 
asymptotic behavior obtained for the evolution of the same wave function for the case of an axially symmetric 
massless scalar in a Minkowskian background 



AG p3/2 1 * ^ 2p2 



2 V 4Cj p-^ 

we see that the leading behavior far from the axis (large p and fixed t) is the same in both cases. 



B. Asymptotic expansions in r 

In principle we only have to consider the situation in which p 7^ because tjj{0, t) — 0. The asymptotic behavior in 
T for the integral 



Jo{pq)Ji{(^q) exp[-jr(l 



|dg 



(23) 



is obtained by the same methods used in previous sections and is 

" 2^^1ogT ''''P + '^^ + +{P + ^) logT])r[-z(p + a)] 

~ exp + cr) - i[i^T + (P + cr) logr])r[i(p + a)] 

+ exp (|(p - <t) + i[-T + (p - cr) logr + |])r[i(CT - p)] 
+ exp {^{<T - p) - i[T + [p - a) log r + T[i{p - a)] } + 0(1/ log' t), 

so that we have that the asymptotic behavior for ^->{p, r) is given by 



/ 2p 
AG{al-al) 



[(72 S'(p, 0-2, t) - (7lS'(p, Cri,T)] 
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This displays the slow decay in time that is characteristic of the system. As before one can compare it with the one 
corresponding to the free, massless, axially symmetric scalar field that is given for large values of r by 



2p 



(p + -2-r)n(^ 



4p(T2 



- {a2 - pf 



4pg-i 

-"-(p-Tl)^ 



- (71 if 
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and decays to zero as 



2t2 



2p 



4G(a| 



This behavior in t means that if the evolution of the initial wave function is such that at some instant of time t 
the probability of finding the particle in the vicinity of the symmetry axis builds up, it will remain high for a large 
interval in r as the asymptotic behavior obtained above shows. In fact this is what happens in this case as can be 
seen in figure Finally it is interesting to notice that the fall-off in the r direction is much faster than the one given 
by the r asymptotic expansion for our system. 

C. Asymptotic expansions for p, a, and t simultaneously large 

We discuss now the obtention of an asymptotic approximation that is valid in a "classical limit" corresponding to 
taking p and t large while keeping their relative values as we did in previous sections. We will also take ci and (72, 
that define the support of the wave function at the initial time, large in comparison with the length scale AG. This 
will allow us to use the same type of asymptotic expansion that we have used in previous sections. The procedure 
should be clear by now so we skip the details here, we just use the contour integral representation introduced above 
for the Bessel functions and write the integrals in (|22|l in terms of them. The asymptotic expansion for H23|l is then 
obtained as the sum of a boundary term contribution that coincides with 



Jo{pq)Jiiaq) exp(-iTg) dq 



plus some extra contributions. The boundary term is 
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and the extra contribution is given by 



e{T-a- p) 



exp 
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^2Tiap{a + p) log ^ 



+ e{T-cr + p)e{cr 
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Introducing these expressions in H22() we finally obtain the desired asymptotic approximation for the wave function in 
the limit when all the lengths arc significantly larger than 4G. 

These asymptotic expansions allow us to explore different possibilities as far as the width and the position of the 
support of the wave function at t ~ to is concerned. We may consider the case in which the support -in the scale 
defined by 4G- is wide or narrow. ^"^ In the first case the wave function evolves in a way that closely resembles the 
propagation of an initial wave function of this type for a free axially symmetric field in a Minkowski background. One 
can easily see that most of the probability amplitude for large values of p is concentrated along the lines t = po + p 
and T — p ~ Po in. the (p, r) plane. These two lines define trajectories that can be interpreted as null geodesies of 
an emergent spacetime metric. Notice that they are defined with a resolution of the order of the width of the initial 
support of the wave function (see figure |SJ). The other case defines a situation when the matter density is in some 
sense high and then displays a behavior that can be interpreted as due to self gravity effects (see figure E)). Also in 
this case, specially when the initial support is close to the axis, there is a build up of the probability amplitude at 
p = that decays subsequently in the very slow fashion characteristic of the model. This means that the probability 
to find the particle in the vicinity of the axis remains high for a long time and, as a consequence, the probability of 
finding it on the "light cone" is much lower. This is shown in figure El 

As we have already discussed in the case of the two-point function and the Newton- Wigner propagator there is 
significant enhancement of the probability to find field quanta close to the axis. This can easily be seen by comparing 
this wave function to the one corresponding to quanta of an axially symmetric massless scalar field in a 2-1-1 dimensional 
Minkowskian background. As the asymptotic behavior of both shows (and the figures clearly display) the probability 
near the axis remains significantly higher in the gravitational case if the support of the initial wave function is narrow. 
Another interesting feature that can be seen is the persistence of a footprint of the wave function in the range of p 
where the support of the initial wave function is. This can be understood by realizing that it is possible to write the 
Newton- Wigner propagator H18|) as the sum of a delta function, multiplied by a time-dependent phase, and a rapidly 
convergent integral. 



In the first part of the paper we have discussed in some detail the exact non-perturbative quantization of Einstein- 
Rosen waves coupled to a massless, axially symmetric scalar field. This provides us with an interesting system where 
gravity is coupled to a matter field that can be exactly quantized and retains some of the features of full general 
relativity. In particular, and even though we have performed a symmetry reduction, we still have and infinite number 
of local degrees of freedom both in the gravitational and matter sector. We also have (some) diffeomorphism invariance 
that has been taken care of through gauge fixing in the Hamiltonian formalism. 

One of the uses of this toy model is to explore the possibility of getting relevant information on the nature of 
quantum space time by using the probe provided by the matter field. This is particularly useful in the setting of 
Einstein-Rosen waves owing to the fact that the object that encodes the gravitational degrees of freedom is a scalar 
field and, hence, its relation with the four dimensional metric -specially in the quantum case- is somewhat indirect. 
Even though we have discussed some specific issues of the combined system (such as two-point functions involving 
both the scalar and gravitational field) we have mostly used the matter field in this role of spacetime probe. 

In order to obtain geometric information on quantum spacetime one must try to work with objects of direct 
physical interpretation in spacetime terms. One such object is the two-point function, in fact, the field commutator 



The Newton- Wigner propagator (-R2|'^{t2, corresponds to the hmit when the support is infinitely narrow. 



VIII. CONCLUSIONS AND COMMENTS 
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FIG. 5: Probability density for a wave function with a wide support at the initial instant of time. The density plot is made 
with the asymptotic approximation discussed in the text to avoid long and slow numerical integrations. Saturated (white) 
areas correspond to divergencies of this asymptotic expansion that are not present in the true wave function. This can be seen 
in the sections plotted in the right hand side of the figure. Notice that in this case the two null curves signalling the light cone 
are well defined -certainly better than in the narrow case shown in figure (S)- and the probability inside the light cone or at the 
axis becomes very low as t grows. 



has already been successfully used in previous work of the authors to understand microcausality 2]. Here we have 
also considered two-point functions with an approximate interpretation (as in ordinary Minkowski space QFT) of 
propagation amplitudes from one spacetime event to another. As is well known this interpretation is only approximate 
because the Hilbert space vectors obtained by acting on the vacuum state with the field operator do not form an 
orthonormal set. In order to overcome this difficulty we have introduced an orthonormal basis of position state vectors, 
labelled by the radial coordinate, and used them to define radial wave functions. These states are a generalization 
of the Newton- Wigner states of ordinary QFT. With the aid of this orthonormal basis we have considered first the 
Newton- Wigner propagator {R2\U{t2, and we have used it later to study the time evolution of a radial position 

wave function. 

The results obtained in all the approaches are different but compatible. They can be summarized as follows: 



• There are some interesting physical effects happening at the symmetry axis. In particular all the approaches 
that we have followed (the approximate ones provided by the two-point functions, the ones given by the Newton- 
Wigner propagator, and the radial wave functions themselves) suggest a significant enhancement of the probabil- 
ity to find field quanta there. In principle one could expect such a behavior due to backscattering -as it happens 
when one considers the classical system-. However we think that the comparison with the massless axially 
symmetric field propagating in a Minkowskian background suggests that most of it is due to a combination of 
quantum and gravitational effects. 
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FIG. 6; Probability density for an initial wave function with a narrow support at the initial instant of time. Notice the enhanced 
probability close to the axis and the position of the light cone structure that appears. 

• The probability amplitudes show an enhanced probabihty of finding scalar field quanta on some lines of slope 1 
in the {p, r) plane that can be interpreted as approximate null geodesies of an emergent metric. Even though we 
cannot obtain other types of geodesies with our massless fields it is reassuring to see the emergence of approximate 
spacetime trajectories with a clear physical interpretation. This approximate null geodesies correspond in this 
case to the ones of the Minkowski metric in 1+1 (or rather 2+1 with axial symmetry). It is worthwhile to note 
that as the one particle states that we are using are "the closest ones" to the vacuum state (other than the 
vacuum state itself) we are seing in an operational way the appearance of a classical fiat spacetime as far as 
(some) of its geometric properties are concerned. 

• Another interesting effect is the persistence of the amplitudes in the support of the initial wave function. This 
is clearly displayed in the behavior of the wave function itself as it evolves in time and shows up in the Newton- 
Wigner propagator. A somewhat similar effect appears in the two-point function. The divergence at pi = p2 for 
these objects is, however, a consequence of the fact that we are not regularizing the fields. Even if we regularize, 
for example by introducing a cut-off 0|, we would observe a large value for this function at pi = p2- 

It would be interesting to apply the methods developed in the paper and use matter field quanta to explore quantized 
geometries for states representing classical configurations corresponding to arbitrary solutions for the Einstein-Rosen 
waves. In particular it would be illuminating to compare the results with those obtained by quantizing the matter 
fields in the curved backgrounds provided by such solutions (where particle creation effects may play a relevant role). 
This would first require us to find suitable semiclassical states for the system describing a non trivial gravitational 
part and a simple matter part and study their quantum evolution. As long as they can be found our approach should 
lead to unambiguous answers to questions related to the emergence of classical trajectories for quantum tests particles. 
We are working on this problem in the present moment. It should be emphasized, however, that a direct comparison 
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between both approaches may be difRcult because of the very different Hilbert spaces used in their quantization. 
In particular the use of Newton- Wigner states in the present scheme is very useful because n-particle subspaces of 
our Hilbert space are stable under the quantum evolution of the system. This allows us to rely on a probabilistic 
interpretation of one-particle wave functions. For a scalar field evolving in a general Einstein- Rosen curved background 
this may be no longer possible due to particle creation effects. 

In our opinion the model provided by Einstein-Rosen waves, free or coupled to cylindrically symmetric matter, is 
an excellent test bed to discuss issues in quantum gravity. Of course there is always the issue as to what extent the 
results obtained are equivalent (at least in a qualitative way) to real effects in a full theory of quantum gravity or 
artifacts of the symmetry reduction. In this respect we think that the effects that we have described above admit a 
sensible interpretation and give interesting hints about the behavior of quantized gravity that we want to explore in 
the future. 
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APPENDIX A: NUMERICAL COMPUTATION OF INTEGRALS INVOLVING TWO BESSEL 

FUNCTIONS 



Throughout the paper we have repeatedly encountered a class of improper integrals involving products of Bessel 
functions. The oscillatory character of the integrand, and the fact that they are defined on the half-real line and 
are usually conditionally convergent makes it necessary to find an efficient way to compute them. We provide such 
a method here. This has been extensively used in the several plots that appear in the paper. The seamless mesh 
between these numerical results and the analytic approximation provided by the asymptotic expansions is an eloquent 
proof of both the accuracy of the asymptotic approximations and the numerical results. 

Let us consider integrals of the type 
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with E Z, and f{q) a sufficiently regular function such that the integral is, at least, conditionally convergent. We 
first change variables according to s = 6""^ to get 
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- Ji^{-P log s)J^{-a log s)/(- log s) ds. 



We now write it as the sum of two contour integrals, h and I2, in the complex s-plane defined on the paths Ci 
{iu : u G [0, 00)}, C2 = {1 + iu : u G [0, 00)}: 



Ii{p,a,T) = e 
h{p,(J,T) = - 



POO 
-IT I 

Jo 

.00 ^ 



u 

-ru 



+ iu 



-J^ii-p log iu)J^{-a log iu)f{- log iu) du , 
J^(-/9log[l + iu])J„{- a log[l + iu])f{- log[l + iu]) du . 



For the functions / that appear in the paper the second integral I2 is very well behaved, because of the exponential 
fall-off of the integrand and its non-singular character. It can be computed numerically without difficulty. On the 
other hand the integrand in 7i has a nasty oscillating behavior in the vicinity of u = although the integral itself is 
convergent. A way to turn it into a much tamer object is to use the contour integral representation for the Bessel 
functions introduced above to write it as the multiple integral 
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where 71,2 are simple closed paths surrounding the origin in C. By choosing these paths appropriately (satisfying, for 
example, the conditions ^e{zi — 1/zi) < and ^e{z2 — I/22) < 0) it is possible in many cases to guarantee that 
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is convergent. This allows us to change the integration order in ljAl|l . Furthermore, for specific choices of the function 
/ (most of the cases appearing in the paper!) this last integral in u can be exactly obtained in analytic form thus 
leaving us with a double integral representation for Ii, defined on a set with the topology of a torus, of a perfectly 
regular function (except for a measure zero set of values of the parameters p, cr, and r). For example if /i = 0, = 1, 
/ = 1 we have 



47r^ 



71 



dzi 
zi 



dz2 

— 2~ 6Xp 



72 



Z2 



1 

Z2 



(*--logr))r 



Z2 



1 

Z2 



The integrand is now well behaved and the integral can be computed numerically in an efficient and quick way by 
ordinary methods. 
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